The effects of k-dependent self-energy 
in the electronic structure of correlated materials 
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It is known from self-energy calculations in the electron gas and sp materials based on the GW ap- 
proximation that atypical quasiparticle renormalization factor {Z factor) is approximately 0.7 — 0.8. 
Band narrowing in electron gas at Ts = 4 due to correlation effects, however, is only approximately 
10%, significantly smaller than the Z factor would suggest. The band narrowing is determined by 
the frequency-dependent self-energy, giving the Z factor, and the momentum-dependent or nonlocal 
self-energy. The results for the electron gas point to a strong cancellation between the effects of 
frequency- and momentum-dependent self-energy. It is often assumed that for systems with a nar- 
row band the self-energy is local. In this work we show that even for narrow-band materials, such 
as SrVOa, the nonlocal self-energy is important. 

PACS numbers: 71.10.-w, 71.15.-m, 71.27.-|-a 
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I. INTRODUCTION 



The electronic structure of a certain class of materials, 
commonly known as strongly correlated materials, is cru- 
cially determined by electron correlations. In these mate- 
rials one-particle theory, predominantly the local density 
approximation (LDA) within density functional theorji, 
is often far from sufficient in providing an accurate and 
reliable description of the electronic structure. The rea- 
son for the failure of one-particle theory may be traced 
back to the presence of many configurations close in en- 
ergy arising from a partially filled narrow band character- 
istic of these materials, similar to the situation in atoms 
with partially filled shell. It is not surprising that de- 
scription of the electronic structure in terms of a single 
Slater determinant is not satisfactory. One of the most 
fruitful approaches in treating the electronic structure 
of these materials is the Green function technique with 
a nonlocal and energy-dependent self-energy. Among 
the Green function approaches, the combination of the 
LDA and the dynamical mean- field theory (DMFT)^, i.e., 
the LDA-f DMFT method^^, is perhaps the most widely 
used technique in describing the electronic structure of 
strongly correlated materials. 

In many calculations based on model Hamiltonians it 
is often assumed that the self-energy is local, meaning 
that it has no momentum or k dependence and the Hub- 
bard U used in solving the impurity problem is assumed 
to be static. How these two approximations affect the 
resulting band structure remains to be investigated. Ex- 
perience with the electron gas teaches us that there is 
a strong cancellation between the effects of k-dependent 
self-energy and frequency-dependent self-energy. This is 
revealed by the fact that effective mass of the electron 
gas with density relevant for most materials is quite close 
to unity. The effective mass in general depends on two 
factors: the derivative of the self-energy with respect to 



frequency, which tends to enhance the effective mass, and 
the derivative of the self-energy with respect with the k- 
vector, which tends to have the opposite effect. In the 
electron gas these two factors tend to cancel each other, 
hence the effective mass being close to unity^'^. We have 
no apriori reason to believe that such an almost com- 
plete cancellation remains true in real materials, espe- 
cially in the strongly correlated systems. Indeed, it is 
often claimed that for narrow-band materials the self- 
energy has little k-dependence so that a local but dy- 
namic theory with frequency-dependent self-energy pro- 
vides a good approximation to the full self-energy. 

The k-dependent self-energy becomes even more rel- 
evant with the recent progress in solving the im- 
purity problem with a frequency-dependent Hubbard 
This new algorithm makes it possible to per- 
form LDA-I-DMFT calculations using a dynamic Hub- 
bard U, instead of a static one as in conventional calcu- 
lations. This raises an important issue concerning the 
role of nonlocal or k-dependent self-energy. The use 
of dynamic U might lead to an overestimation of band 
narrowing or mass enhancement because the dynamic U 
produces a self-energy that has a stronger frequency de- 
pendence or larger derivative at around the Fermi level. 
One might anticipate qualitatively that the use of a dy- 
namic U ought to be counterbalanced by the inclusion 
of a k-dependent self-energy to restore the correct mass 
enhancement. 

The purpose of this paper is to investigate the effects of 
momentum- and frequency-dependence of the self-energy 
on the electronic structure of materials with partially 
filled narrow bands. To address the forementioned is- 
sues, instead of using a model Hamiltonian, we will em- 
ploy the GW approximation (GWA)iSiii which allows us 
to perform accurate calculations on real materials. As a 
test material we study specifically the electronic struc- 
ture of SrVOa which has been widely studied in the lit- 
erature both experimentally by means of photoemission 



2 



spectroscopjsiSrJ^ and theoretically^"—. We believe the 
results are general and will be directly relevant to realis- 
tic electronic structure calculations. 



II. THEORY AND METHOD 

The correlation part of the self-energy (excluding ex- 
change) in the GWA is given b y^°'^^ 



E^(r, r'; u) ^ i f — G(r, r'; u + J)W^{v, r'; u') (1) 



where 



/oo 
-OO 



,r(r,r';c.') 



(9) 



W" -V. 



(2) 



In this work, the one-particle band structure {V'kr! ,ekn} 
is taken to be the LDA one. 

The full self-energy may be expanded in terms of Bloch 
states V'kn: 

S(r,r';u;) = ^ V'k„(r)S„„,(k,c.)^*„(r'). (10) 

knn' 

The Bloch states may be expressed in terms of some Wan- 
nier orbitals {(ySRn} 

■!/'k«(r) =^exp(-zk • R)(^R„(r) (11) 

R 



V is the bare Coulomb interaction and W is the fully 
screened interaction calculated within the random-phase g^^^j ^-^^ self-energy becomes 
approximation (RPA). Using the spectral representation 
of G and W 



G(r, r'; oj) V^kn(r)V'i:„(rO ^ V'k„(r)^;:„(r') 

^ UJ - Ekn -i5 ^ 



^ - Ekri + 



(3) 



S(r,r';aj) 

= ^ ^exp[-ik- (R-R')]¥'R«(r)S„„'(k,w)<^R'„'(r')• 
(12) 



knn' R,R,' 



Vr=(r,r';w) = / 



where 



uj — u)' — i6 



duj 



The local self-energy centered at lattice site R is defined 

„, I ,N according to 
, S(r,r';a;') 



Biuj) = ImM^^(w)sgn(cj), 



w — w' + 

(4) 



(5) 



E^^(r,r';L.) = ^ <^R„(r)S„„- (k, ^)^r„. (r') 

knn' 

= ^^R„(r)S;-(..)^^,„,(r'), (13) 

nn' 

where 



the spectral function of the correlation part of the self- 
energy can be expressed in terms of the imaginary part 
of the screened interaction as follows^: 



(14) 



Since Y}^ is independent of R we may choose R = and 
r(r,r';w < \x) '0kn(r)S(r, r'; ek„-w)V'k„(i"')^(£k«-w),^i.ite 



kn 



(6) 



Ei-(r, r'; 0.) = ^ ^n(r)Sj,-, (..)<, (r') (15) 

occ nn' 

r(r,r';a. > ^ ^ Vk„(r)i?(r, r'; a.-£k„)^L(r')^(^-ek„),^here it is understood that 

kn 

(7) 



where 



<^n =F (fOn- 



(16) 



I The above formulation is quite general and we now fo- 

r(r,r ;uj) = Iml]''(r, r ;i:j)sgn(w - fi). (8) cus on SrVOa which has a narrow band, derived from 

TT .... 

the t2g orbitals of the vanadium, crossing the Fermi 
The real part of the correlation self-energy (excluding level. The t2g band is well separated from the rest of the 



exchange) is given by the Hilbert transform 



bands so that we need only consider Wannier orbitals 
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constructed from Bloch states belonging to this band. 
Following the method of Marzari and Vanderbil t^^'^^ , 
we construct the (maximally localized) Wannier orbitals 
fRn according to 



k 



ik-R 



km/ 



km/ 



(17) 



(18) 



iV'kn) is a one-particle Bloch state which may be chosen 
to be that of the LDA. The matrix is a square matrix 
when the subspace is well separated but for entangled 
bands, S is not necessarily a square matrix, the num- 
ber of band index n may be larger than m. The matrix 
S is chosen to make the Wannier orbitals as localized 
as possible according to a prescription by Marzari and 
Vandcrbilt^i^. 

Using the expression in the expectation value of 
the local self-energy in a Bloch state V'kn is given by 



(^k„|s'°^(t^)Vk„) = 



E 

mm'Ct2g 



(19) 

where is a maximally localized Wannier orbital de- 
fined in (|17p for the central cell R = and the matrix 
elements Sjjim' taken in the Wannier gauge i/j"^^ de- 
fined in Eq. (fTSl) : 



k 



(20) 



is restricted to the orbitals 



The sum in Eq. ([T 
because the Bloch state i/'kn belonging to the t2g band 
has no component outside the t2g subspace. In the case 
of SrVOa, due to cubic symmetry, the self-energy in Eq. 
()20|1 is diagonal and independent of the orbital index m 
so that the matrix element of the local self-energy in the 
Bloch state V'kn in Eq. (fT9|) becomes independent of both 
k and the band index n. 

The quasiparticle band structure is obtained from 



Ekn Ekn + (^knl ReS(£:kn) " Wj;c|'0kn) 
~ S-kn I ^kn 



where 



9ReE„n(k,a;) 
duj 



(21) 



(22) 



is the quasiparticle weight. The angle-resolved spectral 
function is calculated as follows: 



ImE„„(k,cj)| 



„ - ek« - Rel]„„(k,a;)] + [Im I]„„(k, w)] 

(23) 

where we have assumed that the self-energy is diagonal 
in the band index, and the total spectral function is given 

by 



A(a.)=5]A(k,c.). 



(24) 



The band structure calculation is based on the 
full-potential LMTO implementationj^I The exchange- 
correlation functional is the local density approximation 
of the Cepeley- Alder type^^. The GW calculation uses 
mixed basis consisting of products of two atomic orbitals 
and interstitial plane wavesi ^^'^° The 8x8x8 mesh is used 
for Brillouin-zone integration. More technical details are 
found elsewherei^ 



III. RESULTS AND DISCUSSIONS 

A. Local vs nonlocal self-energy 

We first examine the results for the real and imagi- 
nary part of the self-energy at some representative k- 
points r, X, and R shown in FigHJ Im S of the occupied 
states are distinctly different from those of the unoccu- 
pied states. The former have a peak structure below the 
Fermi level at around —3 eV while such structure is es- 
sentially absent for the latter. The structure of Im E can 
be understood by taking the matrix element Im S in Eq. 
(El) and m-. 



^^qm/ 



rm(q,i:^ < M) = ^ ('0qmV'kn|S(£kri " w)|?Akri'0q 
kn 

X eiskn - w), (25) 

unocc 

r,„(q,i:j > ^) = ^ {lljq„l1pUn\B{u - £k«|V'kn'0q 



f^qm/ 



kn 



X 9{UJ - Ekn)- 



(26) 



Since B is proportional to Im W the structure in Im S 
is essentially determined by the structure in Im W with 
intensity governed by the overlap between the state ipqm 
and the occupied or unoccupied states ■(/'kn ■ Whether the 
structure in Im W is carried over to Im E depends on the 
character of the state V'qm ■ If the state ipcim is occupied, 
there will be a strong overlap with occupied states V'kn 
so that the intensity of r(a; < ^) may be expected to 
be stronger than the intensity of r(a; > fi). Indeed, the 
peak structure in Im S for the unoccupied state at the X 
and R point is stronger above the Fermi level than below. 
This structure can be traced back to the peak structure 
in Im W at around 1 — 2 eV as can be seen in Fig [51 
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FIG. 1: (Color online) Self-energy at F, X and R points. The 
local self-energy is also shown. Upper panel shows the real 
part of the self-energy. The straight line represents y = ui — 
Ekn, where tkn is the LDA eigenvalue. Lower panel is the 
imaginary part. 
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FIG. 2: (Color online) Imaginary part of the screened 
Coulomb interaction in the Wannier basis. 



Re E of the occupied states have more structure around 
—3 eV than for the unoccpied states. This structure origi- 
nates from the peak structure in Im S at around the same 
energy. The straight lines represent 



y = W - ffkn 

whose intersection with or proximity to the real part of 
the self-energy below or above the quasiparticle energy 
signals the formation of satellites, provided the imagi- 
nary part of the self-energy at the intersection energy is 
sufficiently small to produce a discernible feature. It is 
clear from the figures that the straight lines only cross 
Re S at one point so that a well-defined satellite feature 
is not expected. However, for the unoccupied states at 
X and R the straight lines come close to the proximity 
of Re S and we expect a formation of a weak satellite at 
around 4 eV above the Fermi level. 

To see the difference between the local self-energy and 
the k-dependent self-energy we plot in Fig. [1] the real 
part of the self-energy at some k-points and compare 
them with the local self-energy, which is the average of 
the self-energy over the Brillouin zone as defined in Eq. 
([20| . The imaginary part of the local self-energy is also 
shown. It resembles the self-energy of the unoccupied 
states because the sum in Eq. ((T4)) is dominated by the 
unoccupied states since the number of unoccupied states 
is much larger than that of the occupied states. The 
contribution from the small number of occupied states 
corresponding to one t2g electron is weighted down by 
the contribution from the rest of the k-points. The real 
part of the local self-energy is shown in FiglT] and similar 
to the imaginary part, it resembles that of the unoccupied 
states with little structure below the Fermi level. A local 
GW self-energy would therefore not produce any distinct 
satellite below the Fermi level at around —3 eV. 



IV. QUASIPARTICLE BAND STRUCTURES 

In Fig. [3] the quasiparticle band structures obtained 
within the LDA, the GWA and the local GWA are com- 
pared. The effect of the k-dependence of the self-energy 
is striking. The band width obtained from the full GW 
calculation is reduced by almost a factor of two when the 
k-dependence is neglected. This result is consistent with 
the following consideration for the effective mass. The 
effective mass to* is given by 



TO 
TO* 



= Z 



1 aRei;(fc,cj) 



de/dk 



dk 



(27) 



where 



Z = 



1 - 



9ReS(fc,cj) 



(28) 
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FIG. 3: (Color online) Quasiparticle bandstructure in LDA, 
GW, and local GW approximation. 



It is clear from the above expression that if the self- 
energy is assumed to be local or onsite then the effec- 
tive mass or the quasiparticle band width is determined 
by the Z factor only since Re E is approximately linear 
within the band width. However, in the case of SrVOa, 
the GW quasiparticle band width is narrowed from the 
LDA band by only 20 % whereas the Z factor is about 
0.5. The discrepancy between the band narrowing and 
the Z factor can be explained by the k-dependence of the 
self-energy. Indeed, when the quasiparticle band struc- 
ture is calculated using a local GW self-energy it is found 
that the band width is reduced by one half from its LDA 
value as can be seen in the figure. Thus, the k-dependent 
self-energy widens the dispersion. This is in accordance 
with the electron gas result, in which the free-electron 
occupied band width is reduced by only 10 % by the 
GW self-energy in the one-shot calculation whereas the 
Z factor is about 0.7. In the self-consistent calculation 
the free-electron occupied band width is widened, rather 
than narrowed. This would be contradictory to the fact 
that the Z factor is ^ 0.7 if the band width were de- 
termined by the Z factor only. Clearly, the k-dependent 
self-energy widens the band width. 

To understand the origin of the band narrowing when 
a local self-energy is used, we plot in Fig. U the differ- 
ence between the full and local self-energies, which rep- 
resents the nonlocal contribution of the self-energy. For 
the occupied states at the F and X points the difference 
is negative whereas for the unoccupied states at the X 
and R points the difference is positive. This implies that 
the occupied states are pushed up whereas the unoccu- 
pied states are pushed down when a local self-energy is 
used, resulting in band narrowing. The effect of band 
narrowing is particularly revealing for the states at the 
X point as can be seen in Fig. U where the nonlocal 




FIG. 4: (Color online) Real part of the nonlocal self-energy 
at F, X and R points. 




FIG. 5: (Color online) k-resolved spectral function at F, X, 
R and M points. 



self-energies for the occupied and unoccupied states have 
different signs. 

Although in the present work we have only presented 
the results for SrVOa we have performed similar calcu- 
lations on other materials such as iron and nickel with 
very similar results. This gives us sufficient confidence to 
believe that the results presented in this paper are quite 
general. 



V. THE SPECTRAL FUNCTIONS 

In Fig. [S]we compare the angle- resolved spectral func- 
tions a.t r, X, R and M points obtained from the LDA, 
the full GW self-energy and the local GW self-energy. 
As can be already expected from the plot of Re S in Fig. 
[T] no strong satellite structure is expected at the F point 
because the straight line y — lo — ekn only crosses Re S at 
one point corresponding to the quasiparticle energy but 
away from this energy it does not come close enough to 
the proximity of ReE. Indeed the (total) spectral func- 
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FIG. 6: (Color online) Total spectral function in LDA, GW 
and local GW approximation. 

tion shown in Fig. |6] only shows a very weak structure 
between —2 and —3 eV with no satellite feature above 
the Fermi level. Experimentally, a satellite at —1.5 eV is 
observed F^^"— 

The spectral functions for the X and R points on the 
other hand show a broad but noticeable satellite feature 
at about 3.5 — 4.0 eV above the Fermi level but with no 
feature below the Fermi level. This is consistent with 
the results for ReE shown in Fig. [T] where the straight 
line corresponding to the unoccupied states come close 
to a peak in Re S around 4 eV. Since there is only one 
d electron, we expect that the satellite above the Fermi 
level corresponding to configuration with two electrons is 
stronger than the one below corresponding to the removal 
of the electron. 

It is known that the GWA tends to overestimate the 
position of the satellite peak arising from a plasmon ex- 
citation. The position of the plasmon peak in W is well 
described by the RPA but the GW self-energy is first- 



order in ly. It is this first-order approximation that 
places the plasmon peak satellite in the spectral function 
too high in energy. It is reasonable to expect a similar 
tendency in low-energy satellites that may be associated 
with the Hubbard bands. Indeed, the experimentally ob- 
served satellite structures have lower energies compared 
with the calculated ones. 



VI. CONCLUSIONS 

The present work reveals that even in materials with 
narrow bands, the nonlocal self-energy has a significant 
effect on the band structure. The band width obtained by 
neglecting the nonlocal self-energy is found to be too nar- 
row compared to the result of a full calculation and con- 
sequently the effective mass is overestimated. Although 
the present investigation has been based on the GWA, it 
is feasible that the result is general. This result indicates 
that calculations based on DMFT using a dynamic U 
where the self-energy is local may overestimate the effec- 
tive mass but further investigations are needed before a 
definite conclusion can be reached because it is not clear 
how the local self-energy in the present work is related 
to the self-energy in DMFT. A promising approach for 
including the nonlocal self-energy on top of the DMFT 
self-energy is the combination of the GW method and 
DMFT'^^ where local correlations responsible for the for- 
mation of Hubbard bands not well described by the GWA 
is taken care of by the DMFT while the nonlocal self- 
energy responsible for modification of the quasiparticle 
energy is accounted for by the GWA. 
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